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Abstract 

Using ReactionKinetics, a Mathematica based package a few dozen detailed mod- 
i-G ' els for combustion of hydrogen, carbon monoxide and methanol are investigated. Es- 

^ \ sential structural characteristics are pulled out, and similarities and differences of the 

O ' 
• i— i , mechanisms are highlighted. These investigations can be used before or parallel with 

^ , 

usual numerical investigations, such as pathway analysis, sensitivity analysis, parame- 
O i 1 ter estimation, or simulation. 
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1. Introduction 



f ■j \ In a previous pap er we presented a Mathematica based program package called 



ReactionKinetics INagy et al 



201211 aimed at symbolic and numerical treatment 



of chemical reactions. The package is especially useful when the number of species 



K^ ' and reaction steps is larger than to allow manual investigations, i.e. if one has dozens 
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or even thousands of species and reaction steps. Since the publication of the previous 
version we made the package capable of reading CHEMKIN files by CHEMKINImport, 
added dozens of new functions such as e.g. 

• CHEMKINExport, 

• MaxFHJWeaklyConnectedComponents, 

• MinFHJWeaklyConnectedComponents, 

• MaxFHJStronglyConnectedComponents, 

• MinFHJStronglyConnectedComponents, 

• FilterReactions, 

and made the package compatible with Version 9 of Mathematica. 

In the present paper we investigate three classes of reactions important both from 
the theoretical and practical points of view in combustion: models of combustion of 
hydrogen, carbon monoxide and methanol, respectively. 

The reaction steps of hydrogen and carbon monoxide combustion form a central 
part of the high temperature combustion of all hydrocarbons and oxygenates. Also, 
hydrogen is an important fuel itself in areas like carbon-free economy, safety issues, 
and rocket propulsion. In the recent years, there has been an increased interest in 
studying the combustion of fuel mixtures consisting of carbon monoxide and hydro- 
gen, referred to as "wet CO" or syngas. These fuels can be produced from coal and 
biomass via gasification, and are considered to be a promising option towards cleaner 
combustion technologies for power generation. Oxygenated organic compounds have 
been proposed as alternative fuels in order to improve the fuel properties and reduce 
particulates and NOx emissions. Methanol is one of the most important oxygenated 
additives since it is the simplest alcohol, has high oxygen content and no C-C bonds. 

The approach we use is absolutely structural in th e sense that n one of results depend 



on the values of the rate coefficients (cf. the paper iBeckl Ill970l0 . We might say that 



we are going to discover possibilities instead of quantitatively dealing with individual 



mechanisms. To put it another way, we are going to raise questions to be answered by 
the chemist, rather than to answer them. 

The structure of our paper is as follows. Section 2 presents the necessary theoretical 
background and describes the models to be investigated. Section 3 shows the results. 

Finally, two electronic supplements are added. First, a Mathematica notebook 
showing all the details of the calculations which may be really useful for those in- 
terested in combustion modeling but annoying for the general audience. Some of the 
resulting figures are also given there. The reader can only redo the calculations if (s)he 



downloads the package itself, from here http : //www . math . bme . hu/~j toth/CES2013 
The date can either be collected from the original authors, or from our database to be 
built in the near future. Second, we also attach the (very long) PDF version of our 
notebook which allows to passively follow what we have done, but this version does 
not need the Mathematica program. (Although the reader can download the program 
CDFPlayer freely, the CDF version of our notebook would not be useful in this case 
because of the heavy use of external data.) 



2. Fundamentals for Formal Kinetics and Combustion 



The basic notions can be found in textbooks such as 



1 1987 



1988]; 



Volpert & Hudyaevl 111 98511 : 



Erdi&T6th[1989]; 



Marin & Yablonski 



1 20 llll etc. 



Feinberg 



2.1. Feinberg-Horn-Jackson graphs, Volpert graphs, Volpert indices of reactions 



Let us consider the reaction 

M 



M 



Y a(m, r)X(m) — > Y f3(m, r)X(m) (r = 1, 2, . . . , R) 

m= 1 m= 1 

with M € N chemical species: X(l), X(2), ..., X(M); R e N reaction steps, 

a(m,r), j3(m,r) e No (m = 1,2, .. .,M;r = 1,2, .. .,K) 
stoichiometric coefficients or molecularities, and suppose its deterministic model 

R 

c„,(t) = f m (c(t)) := ^]06(m, r) - a(m, r))w r (c(t)) 



(1) 



c m (0) = c°eR+ (m=l,2,...,M) 



(2) 
(3) 



— describing the time evolution of the concentration vs. time functions t h-> c m (t) :- 
[X{m)](t) of the species — is based on mass action type kinetics: 

M 

w,.(c) := k r c aU) := k r \\ c a p {p ' r \r = 1,2,...,/?). 

P =\ 

((O is also called the induced kinetic differential equation of the reaction dTJ- — Most 
of the statements does not utilize the special mass action form of the kinetics.) The 
number of complexes N is the number of different complex vectors among a(.,r) and 
/?(., r), i.e. formally it is the cardinality of the set 

{a(.,r);r= 1,2, ...,R] U {J3(., r); r = 1,2,...,/?}. 

The Feinberg-Horn- Jackson graph (or, FHJ-graph, for short) of the reaction is ob- 
tained if one writes down all the complex vectors (or simply the complexes, the formal 
linear combinations on both sides of (Q}) exactly once and connects two complexes 
with a directed edge (or two different edges pointing into opposite directions) if the 
first one is transformed into the second by a reaction step. Let us denote the number of 
weakly connected components (also called linkage classes) of this graph by L. 

The stoichiometric space is the linear subspace of R M generated by the reaction 
vectors: span{y8(.,r) - a(.,r);r = 1,2, ...,/?}; its dimension is denoted by S. Finally, 
the nonnegative integer 6 := N - L - S is the deficiency of the reaction ([T). The 
Fe inberg-Horn-Ja ckson graph will also be used to decide with the method proposed 



by Feinberg 1 1989] if the reaction is detailed balanced or not. 

The Volpert graph of the reaction is a directed bipartite graph, its two vertex sets 
are the species set and the set of reaction steps, and an arrow is drawn from species 
X(m) to the reaction step r if a(m, r) > 0; (species X(m) is needed to the reaction step r 
to take place) and an arrow goes from reaction r to species X(m) if >8(m, r) > (species 
X(m) is produced in the reaction step r). Sometimes it is worth labeling the edges with 
a(m, r) and j3(m, r), respectively. 

It is very useful to assign indices to the vertices of the Volpert graph. This goes 
in the following way. A subset of species is selected, this will be the initial set, (in 
real applications this will be the set of species with positive initial concentrations) and 
the elements of this set receive index zero together with all the reaction steps which 



can proceed once the initial species are present. Next, species without an index which 
can be produced by the indexed reaction steps receive the index one, and reaction steps 
without index which can proceed receive also one, and so on. 

As the Volpert graph is finite, the procedure finishes in a finite number of steps. 
At the end either all the vertices receive an index, or some of them do not receive a 
finite index then one assigns index oo to the vertices without an index. One of the many 
possible interpretations of the meaning of a finite index k is that the given species 
or reaction step can only appear in the Kth step or at the /rth level. In accordance 
with this, species with an infinite index cannot be produced, reaction steps with an 
infinite index cannot proceed with the prescribed initial species of the reaction. This 
st atements and some o t her no t less important ones ca n be found in a precise form e.g. 
in 



Volpert & Hu dvaev 



1 1985] or in the original paper 



Volpert [1972H. A n application 



of the Volpert index in the decomposition of overall reactions is given in 



Kovacs et al. 



1 2004 1. 

Simple examples from combustion theory follow to show the meaning of the defi- 
nitions. 



Example 1 (Mole reaction) The earliest combustion model (which has been given a 
detailed treatment from the point of view of the q ualitative theory of differential equa- 
tions) is probably the Mole reaction ) 







Mold IU936I1 . see also Fig. Q] 
X X + Y — >2X + 2Y. 



(4) 



Example 2 (Robertson reaction) The model proposed in iRobertsonl II 196611 contains 
three species, its FHJ graph is 



A 



B 2B 



B + C — > A + C, 



(5) 



the complexes are A, B, 2B, B + C, A + C, the deficiency isN - L-S =5-2-2= 1. 
The Volpert graph of this reaction is shown in Fig. [2] 

Suppose one takes A as the only initial species, then A and the reaction step 
A — » B gets zero index, B and the reaction steps 2B — > B + C receives 1, finally 



C and the reaction step B + C — > A + C is assigned 2. (Upon selecting B one gets a 
similar result.) However, if one chooses C as the single initial species then all the other 
species and all the reaction steps will have an infinite index. 

We have done all the calculations of the characteristic quantities of reactions us- 
ing the package ReactionKine tics developed in Math ematica and show n also at 
MaCKiE2011 



Napv et al 



2012] and described in detail in 



Toth et al. 



1 201J. Figures 



have also been drawn by the package. 



{Y -> 0, 0) 



{0,1 



|Y,0) 



{X, 2) 




-> Y, 1) 



X, 1| 



•0, 2) 



(X + Y -> 2X + 2Y, 2) 



Figure 1 : The Volpert graph of the Mole reaction (0} 



2.2. Selected models of combustion of hydrogen, carbon monoxide and methanol 

The simplest chain branching combustion reaction, the oxidation of hydrogen is 
already a much more complex system than the Mole and Robertson reactions discussed 
above. Some pe ople feels that this system is well known and no more research is 



needed. However, 



Zselv et al. 



1120 1311 showed recently in a comprehensive mechanism 



comparison paper that the case is not this. The description of the experimental data is 
still not satisfactory and some of the recently published rea ction mechanisms perform 
worse than older ones. Similar comparison was done by 



Olm et al. 



M2013H for the 



oxidation of carbon monoxide. In this work we utilize the mechanism collection of 
these papers, but focus on the structural differences of the mechanisms. Extending the 
investigation with some detailed methanol mechanisms we show that that the suggested 
formal mathematical handling is still applicable for even larger kinetic systems. The 
phenomena are more and more complex as we proceed from hydrogen through carbon 
monoxide to methanol. Correspondingly, the models are larger and more and more 
diverse. 





<C, l| 


[B+C -> A + C. 1} 






{B,0| 


|2 B -> B + C, 0| 


/ 










|2B -.B + C, m| 



Figure 2: The Volpert graphs of the Robertson reaction (0 with B and C as initial 
species, respectively 



3. On the structure of the selected combustion models 

Even the simplest models for combustion usually contain dozens of species of reac- 
tions steps, therefore we can only show some of the results here. E.g. we do not show 
Volpert graphs of the investigated reactions here, because their figures are not useful, 
they can only be used for calculating the Volpert indices. However we have shown the 
Volpert graphs together with Volpert indices above for two simple reactions. 

3.1. Hydrogen 

As a starting point the basic data of the investigated mechanisms are presented in 
Table Q] 

Table 1 : Basic data of the investigated hydrogen combustion models 



Mechanism 


reference 


M 


R 


S=N-L-S 


Ahmed2007 


Ahmed et al. T20071 


8 


38 


29-11-6 = 12 


Burke2012 


Burke etal. [20 121 


8 


38 


31 -12-6 = 13 


CRECK2012 


Healv et al. [20101 


8 


37 


29-11-6 = 12 


Dagaut2003 


Daeaut et al. T20031 


8 


42 


31 -12-6 = 13 


Davis2005 


Davis et al. T20051 


8 


40 


31 -12-6 = 13 


GRI30 


Smith et al. 


8 


40 


31 -12-6 = 13 


Hong2011 


Hone etal. T20111 


8 


40 


31 -12-6 = 13 


Keromnes2013 


Keromnes et al. T20131 


9 


42 


32-12-7 = 13 


Konnov2008 


Konnov T20081 


8 


42 


31 -12-6 = 13 


Li2007 


Li et al. T20071 


8 


38 


31 -12-6 = 13 


NUIG2010 


Healv et al. [20101 


8 


38 


31 -12-6 = 13 


OConaire2004 


6 Connaire et al. [20041 


8 


38 


31 -12-6 = 13 


Rasmussen2008 


Rasmussen et al. [2008a] 


8 


40 


31 -12-6 = 13 


SanDiego2011 


Research T201 11 


8 


42 


31 -12-6 = 13 


S axenaWilliams2006 


Saxena & Williams T20061 


8 


42 


31 -12-6 = 13 


Starik2009 


Starik etal. [20101 


9 


52 


41 -16-7 = 18 


Sun2007 


Sun et al. [20071 


8 


40 


31 -12-6 = 13 


USC2007 


Wane et al. 


8 


40 


31 -12-6 = 13 


Zsely2005 


Zselv et al. T20051 


8 


42 


31 -12-6 = 13 



Now let us start finding the reasons why we have different numbers in different mech- 
anisms. 

Number of species, M All mechanisms contains the same, core set of species: H, H 2 , 
H 2 0, OH, H 2 2 , H0 2 , 2 , O. The Keromnes2013 mechanism, formally, con- 
tains hv as a species, but this is only a notification for the photoexcitation in a 
photochemical reaction. This mechanism is the only one which contains excited 
OH species to describe some ignition delay experiments better. 

There is an other mechanism (Starik2009) which contains an additional species, 
the ozone. It is quite unique to include this species in a reaction mechanism 



intended to be used for the description of combustion processes. 

Number of reaction steps, R The number of reaction steps, R, varies between 37 and 
44, except Starik2009, where this number is 52. Here, (excep t the men tioned 



case) as C. K. Law reported in his comprehensive review paper IILaw , 



200711 the 



number of reactions is approximately 5 times larger than the number of species. 

Deficiency, 6 The number of complexes, N, varies between 29 and 32, except Starik2009, 
where this number is 41 . The number of weakly connected components (or link- 
age classes) is either 11, or 12, except again Starik2009, where this number is 
16. The preliminary data suggest that Starik2009 is structurally larger or richer 
than the other models. 

The deficiencies are large, neither the zero deficiency theorem, nor the one defi- 
ciency theory can be applied. 

Weak reversibility and acyclicity None of the reactions have an acyclic Volpert graph, 
as all the reactions, except CRECK2012, are fully reversible. Accordingly, all 
the reactions are weakly reversible, except again CRECK2012, which has a sin- 
gle irreversible step: H 2 2 + O — > H0 2 + OH. 

3.1.1. Representations of mechanism classes 

From our — let us emphasize: structural — point of view not all the mechanisms in 
Table Q] are different, one has classes with exactly the same structure if the values of 
reaction rate coefficients are disregarded. The classes are shown in Fig. [3] 

During the mechanism development one of the first steps is the decision of which 
species should be included in the mechanism. After this the reaction steps and the 
best possible (as the authors know) set of parameters are selected. This last step, the 
assignment of the rate parameters forms the largest part of a mechanism development 
work. However, we have to keep it in our mind that the parameter set corresponds to 
the previously fixed structure of the model. Therefore, it is important to compare the 
mechanisms from the point of view of their structures. This figure is a demonstrative 
example that the currently published hydrogen combustion mechanisms are different 
not at the level of the parameters, but that of their structures. It is interesting to see 



10 



that a significant number of reaction mechanisms in this collection kept the structure of 
the old GRI1999. It is also interesting, that when the mechanisms are updated most of 
the authors do not modify their structures (see the reaction mechanisms coming from 
the same research group, e.g. SaxenaWilliams2006 and SanDiego2011, or Li2007 
and Burke2012, or OConaire2004 and NUIG2010. The structural relationship of the 
OConaire2004, Li2007, NUIG2010 and Burke2012 is obvious, as they are based on 



some older reaction mechanisms of Dryer's group IMueller et al 



1999]. 



Sun2007 



Davis2005 



OConaire2004 



NUIG2010 





Rasmussen2008 



Iong2011 SaxenaWilliams2006 

o— 



Konnov2008 



Burke2012 



SanDiego2011 



-o 



Keromnes2013 



o 



CRECK2012 



O 



Starik2009 



o 



Dagaut2003 



Ahmed2007 



o 



Figure 3: Classes of hydrogen combustion models. Similar models are connected with 
and edge of the graph. 

This finding forces us to choose a single mechanism from the classes, and we re- 
mind the reader that from now on Davis2005 also represents GRI30, Hong201 1, Ras- 
mussen2008, Sun2007, USC2007; and Burke2012 also represents Li2007, NUIG2010, 
OConaire2004; and Dagaut2003 also represents Konnov2008, Zsely2005; and SanDiego201 1 



11 



also represents SaxenaWilliams2006; whereas each of Ahmed2007, CRECK2012, Keromnes2013 
and Starik2009 form a class alone. The representatives have been selected by the 
caprice of the alphabet. 

3.1.2. Ahmed2007 

As an illustration we show that Feinberg-Horn-Jackson graph in the case of the 
first model. 

We also show the Volpert indices of the species given that the initial species are H 2 
and 2 in Table [2] (The Volpert indices of the reaction steps will only be given in the 
case of Burke2012 below, as an illustration.) 

Table 2: Volpert indices of the species in Ahmed2007 



index 


species 




1 

2 


2 ,H 2 
O, H0 2 , H 
H 2 2 , OH, H 2 



One can also investigate the maximal subgraph of the Feinberg-Horn-Jackson 
graph. It turns out that this substructure is rather stable: The maximal subgraph of 
Ahmed2007, Davis2005 and SanDiego2011 and Starik2009 are the same. Also, in 
the case of Burke2012, CRECK2012 and Keromnes2013 we get the same subgraph. 
And finally, Dagaut2003 (and Konnov2008 and Zsely2005) is special, it is a kind of 
enlargement of the previous graphs. Chemically, OH has four different channels to be 
transformed as opposed to three in the other models. (Maximal can be defined either 
by the number of vertices or by the number of edges, in the special case of combustion 
models selected by us we arrive at the same subgraph using both methods.) 

Let us mention that the three most important radicals in combustion (OH, O and H) 
form a full triangle in most of the maximal graphs except the second series: Burke2012, 
CRECK2012 and Keromnes2013 (including naturally the models represented by these, 
as well.) We shall show below one representative of each type. 
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3.1.3. Burke2012 

The following tables show the indices of species and reaction steps if the initially 
present species are those with zero index. 

Table 3: Volpert indices of the species in Burke2012 



H202 — . _^^ 




Ahmed2007 




^^- H2 + 02 

H20 + H02 
~~ 2H 
H02 + 


~~— 2 OH 


H + H02 


\ 


// 


H20 + 






H202 + O 




H02 + OH 


H20 + 02 






H20 1 




H + OH 


H2 + 


H02 1 




H + 02 


O + OH 










H2 + H02 




H + H202 
2H02 


H20 + OH 


H202 + 02 ^ 




H202 + OH 








~" 20 


02 ^ + 




H2C „ 














H2 + OH 




H + H20 


02 + OH 





Figure 4: The FHJ graph of Ahmed2007 



13 



index 


species 




1 

2 
3 


2 ,H 2 
O, H0 2 , H 
H 2 2 , OH 
H 2 



Table 4: Volpert indices of the reaction steps in Burke2012 



index 


reaction steps 





2 — > 20, H 2 + 2 — > H + H0 2 , H 2 — > 2H 


1 


H + 2 — * + OH, H0 2 +0 — > 2 + OH, H 2 + — ► H + OH, 
H + O — > OH, 20 — > 2 , H + H0 2 — * 20H, 
H + HO z — * H 2 + 2 , 2H0 2 — > H 2 2 + 2 , H + 2 — > H0 2 , 
H0 2 — > H + 2 , H 2 + H0 2 — > H + H 2 2 , 2H — > H 2 


2 


H0 2 + OH — > H 2 + 2 , O + OH — » H + 2 , 2 + OH — > H0 2 + O, 
H0 2 + OH — ► H 2 2 + O, H 2 2 +0 — > H0 2 + OH, 20H — > H 2 + O, 
H + OH — > H 2 + O, OH — > H + O, H 2 2 + OH — > H 2 + H0 2 , 
2 OH — > H + H0 2 , H 2 2 + 2 — > 2 H0 2 , H + H 2 2 — > H 2 + OH, 
H + H 2 2 — * H 2 + H0 2 , 2 OH — > H 2 2 , H 2 2 — > 2 OH, 
H 2 + OH — > H + H 2 0, H + OH — > H 2 


3 


H 2 + 2 — > H0 2 + OH, H 2 + O — > 2 OH, H 2 + H0 2 — > H 2 2 + OH, 
H 2 + OH — > H + H 2 2 , H + H 2 — > H 2 + OH, H 2 — > H + OH 



Let us note that water only appears at the third level. This is the same with some 
models in other classes: CRECK2012 and Kereomnes2013; and in all the other model 
it appears (together with all the other species and reaction steps) earlier, at level 2. 

3.1.4. CRECK2012 

Table 5: Volpert indices of the species in CRECK2012 
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index 


species 




1 

2 
3 


2 ,H 2 
O, H0 2 , H 
H 2 2 , OH 
H 2 



CRECK2012 contains a single irreversible step: H 2 2 + O — > H0 2 + OH. Upon 
going through all the hydrogen combustion mechanisms it turns out that no other model 
contains any irreversible steps. 

3.1.5. Dagaut2003 

Here are the Volpert indices of the species given that the initial species are H 2 and 
2 in Table 

Table 6: Volpert indices of the species in Dagaut2003 



index 


species 





2 ,H 2 


1 


O, H0 2 , H 


2 


OHEX, H 2 2 , OH 


3 


H 2 



The maximal subgraph of the Feinberg-Horn-Jackson graph of Dagaut2003 can be 
seen in Fig. [5] 

3.1.6. Keromnes2013 

Here are the Volpert indices of the species given that the initial species are H 2 and 
O z in Table |7] 

Table 7: Volpert indices of the species in Keromnes2013 



index 


species 




1 

2 


2 ,H 2 
O, H0 2 , H 
H 2 2 , OH, H 2 
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The maximal subgraph of the Feinberg-Horn-Jackson graph of Starik2009 can be 
seen in Fig. [6] 

Another interesting application can be found. Both from the point of view of ther- 
modynamics and from the point of view of reducing the number of reaction rate co- 
efficients one can require that a model be detailed balance (naturally, under the as- 
sumption that temperature and pressure are constant). Applying the pair of conditions 



formulated by 



Feinberg [ 1989] we get the following necessary and sufficient condition 




Figure 5: The maximal subgraph of the Feinberg-Horn-Jackson graph of Dagaut2003 
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in the case of the Keromnes2013 model (ant the models represented by it): 

^26^27^42 = ^25^28^41, fcl£l3&39 = kskukjQt klkukn - k\k\j,k\%, 

k-ik%k\\ — ki,k-jk\2, k 2 k^k$ — k\kik(„ k 2 kt\kc)k 2 \ = fci&3&io&22» 
k 2 hk<)ki9 - kik 7 k w k 2 o, k 2 k 9 k XA k X5 - ^1^10^13^16, 
k 2 k^kx)k\^k 2 ^ky^ — kik-jkiok^k^k^g, k 2 k%k 9 k\ I \k 2 4k-ic — kik-jkioki^k 2 ^k^, 
kik & ki ki3k 24 kl 3 - k 2 k-,k 9 knk 2i k\ A , k x k % kgk X o,k 2A k\ 9 - k^kiokuk^k^, 

9 9 9 9 9 

k 2 k%k<jk\j,k 2 A r k 25 — k\k-ik\Qk\4k 2 T,k 26 , k 2 k 4 kgkgknk 2 4k 3l — kik^kjkiokuk^k 



H202 



H20 + 




17 
Figure 6: The maximal subgraph of the Feinberg-Horn-Jackson graph of 
Keromnes2013 



Further investigations may also used the extended theory bv lGorban & Yablonskvl 11201 ill . 



3.1.7. Starik2009 

Here are the Volpert indices of the species given that the initial species are H 2 and 
2 in Tablel] 

Table 8: Volpert indices of the species in Starik20009 



index 


species 




1 

2 


2 ,H 2 

3 , O, H0 2 , H 

H 2 2 , OH, H 2 



The maximal subgraph of the Feinberg-Horn-Jackson graph of Starik2009 can be 
seen in Fig. [7] 





H20 + 








jw^~~~ 2 ° H - 


H + H02 -_^^ 




H202 ~- 






~~~~~~ H2 + 02 



Figure 7: The maximal subgraph of the Feinberg-Hom-Jackson graph of Starik2009 



3.1.8. Similarities and differences between the models of hydrogen combustion 

We can further analyse the similarities and differences between the class repre- 
sentatives of hydrogen combustion models. One can easily determine which are the 
reaction steps present in one model and missing in the other. This is a very huge table, 
a simpler one is obtained if one only calculates the number of (irreversible, as usual) 
reaction steps present in one model and missing in the other. 

Table 9: Number of different reaction steps in the different models 
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Ahmed 


Burke 


CRECK 


Dagaut 


Davis 


Keromnes 


SanDiego 


Starik 




2007 


2012 


2012 


2003 


2005 


2013 


2011 


2009 


Ahmed2007 





2 


3 








2 








Burke2012 


2 





3 

















CRECK2012 


2 


2 





2 


2 


2 





2 


Dagaut2003 


4 


4 


7 





2 


4 


2 


2 


Davis2005 


2 


2 


5 








2 








Keromnes2013 


6 


4 


7 


4 


4 





4 


4 


SanDiego2011 


4 


4 


5 


2 


2 


4 





2 


Starik2009 


14 


14 


17 


12 


12 


14 


12 






Table [9] shows again that Starik2009 and Keromnes2013 contain a quite few reac- 
tion steps missing in the other models. 

Let us see a single example. The reaction steps what Ahmed2007 contains but 



Burke2012, CRECK2012 and Keromnes2013 do not are H + HO, 



H 2 + O in all 



cases, and in the case of CRECK2012 H 2 2 + O < — H0 2 + OH, as well. Reaction 
steps contained in CRECK2012 and missing in Ahmed2007 are H0 2 ^^ O + OH 
and H 2 2 + O — > H0 2 + OH. The examples should make clear why Tabl^S] is not 
symmetric. 

This shows that in Table |9]reaction steps may mean either an irreversible step, or a 
reversible pair. 

So much should suffice for those interested in the methods, those interested in real 
combustion chemistry should consult all the details in the supplement. 

3.2. Carbon monoxide 

Let us collect now the basic data of carbon monoxide combustion mechanisms into 
Tabled 

Table 10: Basic data of the investigated carbon monoxide combustion models 
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Mechanism 


reference 


M 


R 


6=N-L-S 


Ahmed2007 


Ahmed et al. [20071 


12 


72 


57 - 23 - 9 = 25 


CRECK2012 


Healvetal. T20101 


11 


60 


49 - 19 - 8 = 22 


Dagaut2003 


Dasaut et al. T20031 


12 


68 


52-21 -9 = 22 


Davis2005 


Davis et al. T20051 


11 


60 


47 - 19 - 8 = 20 


GRI30 


Smith et al. 


12 


74 


57 - 23 - 9 = 25 


Keromnes2013 


Keromnes et al. T20131 


12 


64 


52-21-9 = 22 


Li2007 


Li et al. T20071 


12 


78 


61-24-9 = 28 


NUIG2010 


Healvetal. T20101 


12 


78 


61-24-9 = 28 


Rasmussen2008 


Rasmussen et al. [2008a] 


13 


88 


66 - 26 - 10 = 30 


SanDiego2011 


Research [20111 


12 


74 


57 - 23 - 9 = 25 


S axenaWilliams2006 


Saxena & Williams T20061 


11 


60 


45-18-8 = 19 


Starik2009 


Stariketal. T20101 


13 


88 


70-28-10 = 32 


Sun2007 


Sun et al. T20071 


12 


66 


52-21 -9 = 22 


USC2007 


Wane et al. 


12 


74 


57 - 23 - 9 = 25 


Zsely2005 


Zselv et al. T20051 


11 


62 


47 - 19 - 8 = 20 



Since the reaction mechanism of the hydrogen oxidation is a submechanism of that 
of the carbon monoxide oxidation all CO mechanisms can be used for the description 
of the combustion of hydrogen. This means that all mechanisms in Table 8 appear in 
Table 1, but in the previous chapter we focused on the hydrogen submechanisms only. 
Although the references show a large overlap with those in Table [T] here we focus on 
the submechanism describing carbon monoxide combustion. Now let us start finding 
the reasons why we have different numbers in different mechanisms. 

Species and their number, M Here one sees a much more diversified picture. The 
core species present in all the mechanisms are CO, C0 2 , H, H 2 , H 2 0, HCO, 
H0 2 , H 2 2 , O, 2 , OH. Species outsides this set can be seen in Table fTTI 

Table 1 1 : Species present only in a given mechanism 
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Particular species 


Mechanisms 


CH 2 


Ahmed2007, Dagaut2003, GRI30, Li2007, 
NUIG2010, Rasmussen2008, SanDiego201 1 
Starik2009, Sun2007, USC2007 


OHEX 


Keromnes2013 


HOCO 


Rasmussen2008 


o 3 


Starik2009 



Number of reaction steps, R The number of reaction steps, R, varies between 6 and 



. The law of C. K. Law reported in his comprehensive review paper OLaw , 



2007] that the number of reactions is approximately 5 times larger than the num- 



ber of species is fulfilled again. 

Deficiency, 6 The number of complexes, N, varies between 45 and 66. The number 
of weakly connected components (or linkage classes) changes between 18 and 
28. The deficiencies are large, neither the zero deficiency theorem, nor the one 
deficiency theory can be applied. 

Weak reversibility and acyclicity None of the reactions have an acyclic Volpert graph, 
as all the reactions, except CRECK2012 and Keromnes2013, are fully reversible. 
The mentioned two mechanism both contain two irreversible steps: 



HCO + HO, 



C0 2 + H + OH and H 2 2 + O 



H0 2 + OH (CRECK2012); 



2 HCO — > 2 CO + H 2 and HCO + H0 2 — > C0 2 + H + OH (Keromnes2013). 

Although, irreversible steps are acceptable modelling tools, there are at least two 
problems with them. If only one way of the reaction is used the negligibility 
of the reverse reaction step may depend on the circumstances and it is possible 
that the mechanism will be used in such conditions where this simplification 
assumption will not be valid. In case both directions are present in a mechanism 
both of their values should change according to the thermodynamic equilibrium 
if they are re-parametrized. 



21 



3.2.1. Representations of mechanism classes 

The situation is much simpler here, in the case of carbon monoxide combustion 
models. From the structural point of view USC2007 is identical to GRI30, NUIG2010 
is identical to Li2007, and all the other mechanism are different. Therefore we do not 
introduce classes of mechanism here. 

3.2.2. Similarities and differences between the models of carbon monoxide combustion 

Table 12: Species present in one CO model and missing in others 1 



Mechanism 


Ahmed 
2007 


CRECK 
2012 


Dagaut 
2003 


Davis 
2005 


GRI30 


Keromnes 
2013 


Li 
2007 


Ahmed2007 





19 


12 


14 


2 


16 


2 


CRECK2012 


7 





6 


7 


7 


6 


6 


Dagaut2003 


8 


14 





8 


8 


8 


4 


Davis2005 


2 


7 











2 


2 


GRI30 


4 


21 


14 


14 





16 


2 


Keromnes2013 


8 


10 


4 


6 


6 





4 


Li2007 


8 


24 


14 


20 


6 


18 





Rasmussen2008 


16 


34 


24 


28 


14 


29 


12 


SanDiego2011 


4 


19 


14 


14 


2 


16 


4 


SaxenaWilliams 
2006 


4 


7 


2 


2 


2 


4 


4 


Starik2009 


16 


35 


24 


28 


16 


29 


14 


Sun2007 


6 


12 





6 


6 


6 


2 


Zsely2005 


4 


9 





2 


2 


4 


4 



Table 13: Species present in one CO model and missing in others 2 
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Mechanism 


Rasmussen 
2008 


SanDiego 

2011 


Saxena 

Williams 

2006 


Starik 
2009 


Sun 
2007 


Zsely 
2005 


Ahmed2007 





2 


16 





12 


14 


CRECK2012 


6 


5 


7 


7 


6 


7 


Dagaut2003 


4 


8 


10 


4 


2 


6 


Davis2005 








2 











GRI30 





2 


16 


2 


14 


14 


Keromnes2013 


5 


6 


8 


5 


4 


6 


Li2007 


2 


8 


22 


4 


14 


20 


Rasmussen2008 





16 


30 


14 


24 


28 


SanDiego2011 


2 





14 


2 


14 


14 


SaxenaWilliams 
2006 


2 








2 


2 


2 


Starik2009 


14 


16 


30 





24 


28 


Sun2007 


2 


6 


8 


2 





6 


Zsely2005 


2 


2 


4 


2 


2 






As an illustration let us calculate the Volpert indices of Zsely 2005 under the as- 
sumption that the species 2 , H 2 and CO are initially present. 

Table 14: Volpert indices of the species in Zsely2005 



index 


species 




1 

2 


2 , H 2 , CO 

O, H0 2 , H, OH, HCO, C0 2 

H 2 Q 2 , H 2 



The maximal components of the Feinberg-Horn-Jackson graphs are the same as 
those found in the case of hydrogen combustion models. The reason of this is, that 
there are not enough carbon containing species in the models to form larger subgraphs, 
which is not the case with methanol models. 
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3.3. Methanol 

Let us start again with the basic data, see Table [T31 

Table 15: Basic data of the investigated methanol combustion models 



Mechanism 


reference 


M 


R 


6=N-L-S 


Aranda2013 


Arandaetal. T20131 


76 


1063 


661 - 187-71 =403 


Klippenstein201 1 


Klippenstein et al. [201 1] 


18 


172 


122 - 42 - 15 = 65 


Li2007 


Li et al. T20071 


18 


170 


121 -42-15 =64 


Rasmussen2008 


Rasmussen et al. [2008b] 


28 


320 


222 - 75 - 24 = 123 


ZabettaHuba2008 


Zabetta & Huoa T20081 


58 


724 


500 - 163 - 54 = 283 



The analysis of these models is much harder. 

Species, classes of mechanisms The number of reaction steps is around ten times that 
of the species here. 

There is a striking similarity of Klippenstein201 1 and Li2007 at the level of num- 
bers. Really, they use the same set of species, and the only difference between 
their reaction steps is that Klippenstein201 1 contains also the step 



CH 3 Q + H 2 2 === CH3OH + HO z 



beyond the common steps. It is in accordance with the statement of the authors 
that they only made a small change on the structure of Li2007. 

Otherwise, methanol models are so different that the question of classes and their 
representation does not even come up. 

Number of reaction steps, R The number of reaction steps, R, ranges between 18 and 

76. 

Weak reversibility and acyclicity None of the reactions have an acyclic Volpert graph, 
as Li2007, Rasmussen2008 and ZabettaHuba2008 are fully reversible, and most 
of the reaction steps of the two other reactions are reversible. The exceptions are 
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that Aranda2013 contains the irreversible steps 

CH 2 OOH — > CH 2 + OH, 

C 2 H 4 + HOCH 2 CH 2 00 — > CH 2 + CH 2 OH + CH 3 CHO, 

CH 2 + HOCH 2 CH 2 00 — > CH 2 OH + CH 2 OOH + HCO, 

H0 2 + HOCH 2 CH 2 00 — > CH 2 OH + CH 2 OOH + 2 , 

CH 2 CHOH + 2 — > CH 2 + HCO + OH; 

whereas Rasmussen2008 contains the irreversible steps 

N0 2 — >N0 2 -,2N0 2 - — >2NO + 2 , 

and all the other steps in all the other mechanisms are reversible. 

Deficiency, 6 There is no room for the application of deficiency either. As even the 
smallest FHJ graph is too large to be shown here, still we shall deal with the 
largest subgraphs of the FHJ graph of the individual mechanisms. Volpert graphs 
will only be used for indexing, and will show some interesting relationships. 

3.3.1. Aranda2013 

The maximal weakly connected components or linkage classes of the FHJ graph of 
Aranda2013 are shown in Fig. [8] 

Starting from CH 9 as initial species all the reaction steps can finally take place 
and all the species will be produced, except those steps where compounds of nitrogen 
occur in the reactant complex. If one takes {CH 2 0, N0 9 , NH 3 } as the initial set, then 
all the steps are capable of taking place and all the species will be produced and the 
largest index is now 4. Are there any species containing two carbon atoms (or C-C 
bonds, as these expressions are synonymous in this case)? Yes, there are more than 
twenty. And — beyond Aranda2013 — it is only ZabettaHuba2008 where species with 
not less than two carbon atoms can be found, see below. 

3.3.2. Klippenstein201 1 

The maximal weakly connected components or linkage classes of the FHJ graph of 
Klippenstein201 1 or Li2007 are shown in Fig. [9] 
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CH20 + CH3 
CH2CH20H y ,,, -, 

J/mk CH3CH20 

\ a ™ + o i/ 

3 + H ^Z^ C2H4 + OH -~-_J^V (/ 



C2H3 + H20 



» CH2CH 

I 



CH3CHOH 
CH2CHO + H2 



CHCHOH + H2 



CH30 + C0 

CH3CO + 

+ C02 



/r _^w CHL1 + C 



H CH2CHOO ^ C2H3 + ° 2 



H CH20H + CO 



CYCOOP * ^^- ^ ■»!< , CH2CO + OH 

-"" —> CH20 + HCO C2H2 + H02 rj CH2CHO + * ' \\ 

H + OCHCHO 1 — ^ ^^\ 

CHCHOH +0 „-„ 



HCCOH + OH 



Figure 8: The maximal weakly connected components of Aranda2013 



Starting from CH 2 as initial species all the reaction steps can finally take place 
and all the species will be produced. The species H 2 2 and all the reaction steps where 
this species is a reactant species (and only those) will appear latest, only at level four. 

3.3.3. U2007 

No wonder that the situation is the same as with Klippenstein201 1. 
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Figure 9: The maximal weakly connected components of Klippenstein201 1 and Li2007 



3.3.4. Rasmussen2008 

The maximal weakly connected component or linkage class of the FHJ graph of 
Rassmussen2008 is shown in Fig. \W\ 

Starting from CH 2 as initial species all the reaction steps can finally take place 
and all the species will be produced, except those steps where compounds of nitrogen 
occur in the reactant complex. If one also adds NO initially, then all the steps are 
capable of taking place and all the species will be produced. If the set of initial species 
is {CH z O, N0 2 ), then the situation is even better: the largest index is now only 3. 

3.3.5. ZabettaHuba2008 

The maximal weakly connected component or linkage class of the FHJ graph of 
ZabettaHuba2008 is shown in Fig. Q~T] 

Again, if the initial species is only CH 2 0, then no nitrogen compounds (and radicals 
etc.) are produced. However, adding NH 3 initially, all the species and reaction steps 
receive a finite Volpert index. 

Are there any species containing two carbon atoms? Yes, there are: CH 2 CO, 
HCCO. 
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3.3.6. Similarities and differences between the models of methanol combustion 

The mechanisms are too large (especially Aranda2013) to present here all the de- 
tails. However, it is possible to show the number of reaction steps present in the differ- 
ent models and missing in the others, see Table [T6l 





Figure 10: The maximal weakly connected component of Rasmus sen2008 




^L^~ 



CN + H20 



Figure 1 1 : The maximal weakly connected component of ZabettaHuba2008 
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Table 16: The number of different reaction steps in the different models 





Aranda 


Klippenstein 


Li 


Rasmussen 


ZabettaHuba 




2013 


2011 


2007 


2008 


2008 


Aranda2013 





897 


899 


759 


737 


Klippenstein2011 


6 





2 


8 


18 


Li2007 


6 








6 


16 


Rasmussen2008 


16 


156 


156 





120 


ZabettaHuba2008 


398 


570 


570 


524 






Thus, 897 in the first row, second column in Table [TBI means that there are altogether 
897 reaction steps enumerated in Aranda2013 but missing in Klippenstein2011. Note 
that the table is not symmetric, it should not be in general. 

Let us see a few examples in more detail. As Li2007 is a proper subset of Klippen- 
stein201 1, there is no reaction step present in the first one and missing in the second. 

The reaction steps present in Klippenstein2011 and missing in ZabettaHuba2008 
are: 



2CH 2 


^=^ CH 2 OH 


+ HCO, 


CH 2 + CO ^^ 2 HCO, 




CH 3 + 


CO ^ CH 3 


+ co 2 , 


CH3O + H ^^ CH 3 + OH, 




2 CO + H 


^2 HCO, 


CH 2 C 


» + H 2 2 ^^ CH 2 OH + H0 2 




CH 2 + H 2 2 


^ CH3O + 


H0 2 , CH3O + H 2 2 ^^ CH3OH + 


H0 2 , 




HCO + H0 2 ^^ 


C0 2 + H + OH 





Finally, Rasmussen2008 contains a few reaction steps among nitrogen compounds (in- 
cluding two irreversible steps) which are not present in the huge Aranda2013. These 
are as follows. 





CH 2 + H 2 2 ^ 


CH3O 


+ H0 2 


HN0 2 ?= 




HO : 


+ N0 2 ^= 


±N0 3 


+ OH, 


NO3 ;=^ N( 




H + NO3 


^N0 2 + 


OH, 


H0 2 + 


N0 3 ^^ NO 


NO 


3 + o^=± 


N0 2 + 2 , 


N0 2 


— >N0 2 -, 2N0 2 



H + NO z , 



NO + 2 , 



N0 2 + Q 2 + OH, 



2NO + 0, 
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4. Discussion and outlook 

The major application of the methods outline in the paper might be that before start- 
ing quantitative analysis like estimating the reaction rate constants etc. one analyses 
the structure of the models at hand. 

A systematic use of Volpert indexing may also serve for selecting a minimal initial 
set os species: the least number of species which is enough for all the reaction steps in 
a given model to occur and for all the species to be produced. 

Another possible application is that one starts from a big model and deletes reac- 
tion routes obeying some restrictions. E.g. one starts from a CO combustion model 
and deletes reaction steps containing C, thus we should arrive at a hydrogen combus- 
tion model etc. The results are only useful if they are exported to a CHEMKIN file, 
CHEMKINExport will serve for this purpose. 

We also hope that our method can and will be applied in metabolism and atmo- 
spheric chemistry research, as well. 
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